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Abstract. We study the weak solvability of a system of coupled Allen-Cahn-like equations resembling 
cross-diffusion which is arising as a model for the consolidation of saturated porous media. Besides using 
energy like estimates, we cast the special structure of the system in the framework of the Leray-Schauder 
fixed point principle and ensure this way the local existence of strong solutions to a regularised version 
of our system. Furthermore, weak convergence techniques ensure the existence of weak solutions to the 
original consolidation problem. The uniqueness of global-in-time solutions is guaranteed in a particular 
case. Moreover, we use a finite difference scheme to show the negativity of the vector of solutions.Weak 
solutions; cross-diffusion system; energy method; Leray-Schauder fixed point theorem; finite differences; 
consolidation of porous media 


1. Introduction 


Porous solids with fluids moving inside are very important to numerous engineering applications in¬ 
cluding the classical soil compaction and consolidation problem in civil engineering and poromechanics, or 
the biomechanics of bones and tissues, consolidation and subsidence control in environmental engineering, 
seepage of polluted liquids leaking from dangerous reservoirs, oil extraction plants and geothermal reser¬ 
voirs; see for instance chapter 6 in [|3J for basic theoretical accounts and [|4j[l7j, |9j] and references cited 
therein for more modern applications. 

A typical unwanted phenomenon in the consolidation context is the occurrence of phase separation 
between fluid-rich and fluid-poor regions in porous media. Indeed, in such a case the porous medium, even 
in presence of an external pressure, could possibly have in its interior dangerous fluid bubbles (8j. 

In this paper, we study a time-dependent Allen-Cahn-like system modelling the evolution of the macro¬ 
scopic strain and fluid density in a porous media which is able to produce steady states exhibiting a strong 
phase separation between fluid-rich and fluid-poor regions; for details see (5-{T0|. The system we are 
studying (referred to as problem (P) in Section 2.1 1 has two mathematically challenging components: (i) 
a coupled flux (a linear combination of strain and fluid density gradients) resembling this way with cross¬ 
diffusion problems (see [211, e.g.) or with thermo-diffusion problems (see Q e.g.); (ii) the polynomial 
structure of the production term. 

Trusting the working techniques from [5J, we apply a variant of the Leray-Schauder fixed point theorem 
to prove the existence of strong solutions to a regularized consolidation problem (see Section ED and 
then employ weak convergence methods for this auxiliary problem to obtain in the limit of the vanishing 
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regularisation parameter local-in-time weak solutions of the original consolidation problem. Under some 
additional restrictions on the model parameters, we show that the weak solutions exist globally in time 
and are negative. We conclude the paper with numerical illustrations of the solution to our problem and 
point out their non-uniqueness at stationarity for critical parameter regimes. We also briefly discuss a few 
mathematical aspects still open in this context. 


2. Problem and results 

In this Section, we introduce the problem we are interested in and state our main results. In Section [2~4] 
we shall discuss the our main physical motivations coming from the porous media physics. 


2.7. Strong formulation of the problem 

If e denotes the strain and m the fluid density of our porous media (say £2) during a given observation 
time interval (say S), then the strong formulation of the problem we are going to study reads as follows: 


d £ * 

— +div(—kiVe — koS/m) = f\(m,e) 


dm 


+ div(—A^Ve — kfWm) = fi{m,£) 


dt 

e(x,0) = £o(x) 
m(x, 0 ) = mo(x) 

e(h,t) = e D (t) 

m(h,t ) = m D (t) 
ds 

f(h,t) = 0 

ox 

dm., 

& ( ' 2 ’ ,)=0 


in£2 x 5, 

in£2 x S, 

in £ 2 , 
in £ 2 , 
inS, 
inS, 

inS, 

inS. 


( 2 . 1 ) 

( 2 . 2 ) 

(2.3) 

(2.4) 

(2.5) 

( 2 . 6 ) 

(2.7) 

( 2 . 8 ) 


We refer to (2.1 >—(2.5 I as problem (P). 

This paper targets at the weak solvability of problem (P). Before stating our main results, we collect the 
assumptions imposed on the data and parameters involved in the model equations. 

H \: The boundary functions Eff t). mj)(t) are negative continuous for all t E S with d t £n\. \d t mo\ < C for 
a positive constant C. 

Hi'. £o, mo E C(£2) with £o < 0, mo < 0. 

7 / 3 : Let M\ , Mi E M sufficiently large. We take 

/1 (r,s), if |r| <Miand| 5 | <Mi 


fi{r,s) := 


0 , 


otherwise, 


(2.9) 


h(r,s) := 


fi(r,s), if |r| <Miand|s| <Mi 


where f\ , fi : M x ! 


0 , otherwise, 

l. For the setting of interest, we consider f \, fi defined by 


( 2 . 10 ) 




HI A Zs 

k l= 0k 2 =0 


® * h m k \ 


2 E R, nf E N,*,• E {0,.. ij = 1,2, 


( 2 . 11 ) 
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namely we take f \, / 2 to be generic polynomials. 

H 4 : k]. I <2 - ^3 and y are strictly positive constants. 

7 / 5 : k 2 < min {k\ As}. 

H 6 : e 0 ,m 0 GC v (ft), v > 0. 

Assumptions H\ , Hi reflect the properties of suitably rescaled and translated mechanical strain and fluid 
density; Hj, H4 are made so that this scenario fits to the setting described in |9j, while H*,, Hq are technical 
assumptions. 


2.2. Notation 

For a function g = g(x,t), d x g (orVg), d t g indicate the partial derivatives with respect to spatial variable 
x and temporal variable t. Let T, l\ ,/ 2 > 0 be fixed values. Define ft := (h,h), S := (0,T], and L := |ft| = 
h~h- For 1 < p < 00 we denote by ZT(ft) the usual Lebesgue space equipped with the norm || • ||i#>(o)- For 
I < p < 00 and k a positive integer, let W k,p (D.) be the usual Sobolev space with the norm || ■ || w *j,( Q ). We 
write H s (Cl) and || • || h s (q) instead of W S,2 (Q.) and |j • Let T D := { l\ }. We denote by V the space 

V := {q> € H l (ft), (p = 0 on 

and recall the equivalence 


II' lliF(n) ~ II • ||v(np (2-12) 

For p £ [l,oo) we denote by L P (S,B ) the usual Bochner space equipped with the norm 

II ‘ II Lp(5;B)j 

for any arbitrary Banach space B equipped with norm || • ||g. 


Let Bo, B 1 be two Banach spaces. Define the space 

W := {v; v G L 2 (S;B 0 ), d t v G L 2 {S',Bi)} , 


(2.13) 


and take, as a particular case, Bq = V, B\ = L 2 (ft). By [22], Proposition 23.23 (ii) p. 422, we have that 

W ^^C(S;L 2 (ft)). (2.14) 

The following compactness result due to [2| will be useful in our context: 

Theorem 1 (Aubin), Let Bq, B B\ be three Banach spaces where Bq, B\ are reflexive. Suppose that Bq is 
continuously imbedded into B, which is also continuously imbedded into B\, and, morover, the imbedding 
from Bq into B is compact. Let W be defined as in ( 2.1 3[ ). Then the imbedding from W into L 2 (S: B) is 
compact. 

Define the space 

F 2,1 (ft x S) := {(p G L 2 (Q x 5), (p t ,(p xi (p xx G L 2 (ft x 5)} , 
then the following imbedding is a consequence of Theorem [T] 

V 2 '\QxS)^^L 2 (S-H l (Q)). (2.15) 


2.3. Main results 
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Definition 2. The couple 

(e,m) G [£ Z , + L 2 (5;y)n// 1 (S;L 2 (Q))] x [m D + L 2 (S;V)nH\S;L 2 (Ll))] 
is called weak solution to problem (P) if and only if the following identities 

(d,e, (p) L 2 (fl) + h (Ve, v<p) i2(n) = -k 2 (Vm, V<p) L 2 (n) + (/i, <p) L 2 (a) , (2.16) 

(d t m, \lf)i2(n.) +k 3 (Vm), V^l 2 ^) = — fo(Ve, Ny)l 2 {sS) + (A; V)l 2 (q.)- (2-17) 

hold for all (<p, t/r) G V x V and for all t £ S. 

Theorem 3 (Existence). Under the assumptions H \ —H\ there exist at least a weak solution to problem (P) 
in the sense of Definition^ 

Theorem 4 (Uniqueness). Assume Tp-Tp a/id 7/ 5 to hold. Then, for any fixed T G (0°°), it exist at most a 
solution to (P) 7>7 the sense of Definition^ 

Theorem 5 (Boundedness and negativity of e and m). Assume Tp-Tp and /p to hold true together with 
A | —A 3 ( cf Section [ 5 ]). Moreover, assume that the functions f\ . fi are negative. Then the solution ( £,777) to 
Problem (P) ;'.v bounded and negative. 

Remark 6. It is worth noting that Theorem [T| Theroem [4] and Theorem [5] were obtained for the case of 
Dirchlet—Neumann boundary conditions. Note that with minimal modifications of the proofs, we can handle 
other kinds of physically relevant boundary conditions (e.g the periodic case or the Dirichlet-Dirichlet or 
the Neumann—Neumann boundary conditions). 


2.4. Application to the consolidation of porous media 


The problem (P) introduced in Section 2.1 as already announced in the introduction, has a relevant 
application to the theory of Porous Media. In this section we give a very brief account of this theory and we 
refer the interested reader to the paper fTOj for a detailed derivation. 

We introduce the one dimensional poromechanical model (see 1101) whose geometrically linearized ver¬ 
sion is connected to problem (P). Kinematics will be briefly resumed starting from the general statement of 
the model cf. 0 The equations governing the behavior of the porous system are then deduced prescribing 
the conservative part of the constitutive law through a suitable potential energy density <f> and the dissipative 
contributions through purely Stokes term. 

Let fi s := [l\ ,If\ C M, with t \ Jf G M, and B\ := M be the reference configurations for the solid and 
fluid components; see ill]. The solid placement % s : B s x M —> M is a C 2 function such that the map 
Xs(-f), associating to each X s G B s the position occupied at time t by the particle labeled by X s in the 
reference configuration If, is a C 2 -diffeomorphism. The fluid placement map : If x I G I is defined 
analogously. The current configuration B t := x s (B s ,t) at time t is the set of positions of the supeiposed solid 
and fluid particles. Consider the C 2 function (j) : B s x M —>■ Bf such that <j)(X s ,t) is the fluid particle that at 
time t occupies the same position of the solid particle 2f s ; assume, also, that </)(•,7) is a C 2 -diffeomorphism 
mapping univocally a solid particle into a fluid one. The three fields % s , fj. and (j) are not at all independent; 
indeed, by definition, we immediately have that X\(f)(X^.t),f) = ^ s (X s ,f) for any X s G B s and 7 G M. 

The Lagrangian velocities are two maps associating with each time and each point in the solid and fluid 
reference space the velocities of the corresponding solid and fluid particles at the specified time. More pre¬ 
cisely, the Lagrangian velocities are the two maps u a : B a x M —> M defined by setting u a (X a .t) := dXa/dt 
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for any X a S B a , where a = s,f. We also consider the Eulerian velocities v a '■ B t x M —> M associating with 
each point r£8, and for each time t E M the velocities of the solid and fluid particle occupying the place x 
at time f; more precisely we set v a (x,t) := u a {Xa l ( x d)d)- 

In studying the dynamics of the porous system, we can arbitrarily choose two among the three fields Xs, 
Xu and 0. Since the reference configuration B s of the solid component is known a priori, a good choice 
appears to be that of expressing all the dynamical observables in terms of the fields x& and 0 which are 
defined on Z? s . 

It is natural to assume that, if the system is acted upon only by conservative forces, its dynamics is 
described by a Lagrangian density L, relative to the solid reference configuration space volume, depending 
on the space variable X s and on time through (in principle) Xs, 0, Xs’ > Xs’ 0^ Xs, and 0. The Lagrangian 
density is equal to the kinetic energy density minus the overall potential energy density accounting for both 
the internal and the external conservative forces. 

Suppose the fluid component of the system is acted upon by dissipative forces. We consider the inde¬ 
pendent variations §Xs and 8 (j) of the two fields Xs and 0 and denote by 8 W the corresponding elementary 
virtual work made by the dissipative forces acting on the fluid component. The possible motions of the 
system, see for instance [7, Chapter 5], in an interval of time (ti,t 2 ) C M are those such that the fields Xs and 
0 satisfies the variational principle 

sTdtf dA s L(^(X s ,l),...,0(X s ,O) = - f 2 8Wdt (2.18) 

Jt\ JB, J t\ 

namely, the variation of the the action integral in correspondence of a possible motion is equal to the integral 
over time of minus the virtual work of the dissipative forces corresponding to the considered variation of the 
fields. 

The way in which dissipation has to be introduced in saturated porous media models is still under debate. 
In particular, according to the effectiveness of the hypothesis of separation of scales, between the local and 
macroscopic level, Darcy’s or Stokes’ effects are accounted for. We refer the interested reader to fl0| for a 
detailed discussion of this issue. In this paper, we consider the so-called Stokes’ effect, i.e., the dissipation 
due to forces controlled by the second derivative of the velocity of the fluid component measured with 
respect to the solid. A natural expression ffO| is 

8 W:=-[ S[vf(x,t)-v s (x,t)]'[5^f(^ f _1 (x,t),t)-5^ s (^“ 1 (x,t),t)] / (k (2.19) 

Jb, 


where 8 xt is the variation of the field Xt induced by the independent variations 8 Xs and 8 (j), and S > 0. 

In order to write explicitly the variation of the action one has to specify the form of the Lagrangian 
density. In the sequel we shall not consider the inertial effects, so that, the Lagrangian density will be the 
opposite of the potential energy d> density associated to both the internal and external conservative forces. 
It is reasonable to assume that the potential energy density depends on the space and time variable only via 


two physically relevant functions: the strain of the solid and a properly normalized fluid mass density [ 101, 


i.e., 


e(X s ,t) := [(Zs(X s ,t)) 2 - l]/2 and m f (X s ,t) := p 0 M(X s ,t))^ (X s ,t) 


( 2 . 20 ) 


where po.f: B( —> R is a fluid reference density. In other words, we assume that the potential energy density 
is a function of the fields ntf and e and on their space derivative m' ( and e'. 

By a standard variational computation, see | [T0| equation (24)], one gets the equation of motion. In this 
framework, we are interested in the geometrically linearized version of such equations: we assume poy to 
be constant and introduce the displacement fields u(X s ,t) and w(X s ,t) by setting 


X s (X s ,t)=X s + u(X s ,t) and 0(X s ,t) = X s + w(X Sl t) 


(2.21) 
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for any X s £ fi s and t £ M. We then assume that u and w are small, together with their space and time 
derivatives, and write 

rn = po,f(l +w'), m := /?7f-p 0 ,f = PoW, £ ~ u', (2.22) 

where sa means that all the terms of order larger than one have been neglected. We then write the equations 
of motion up to the first order in u, w, and derivatives: 


de 


vde 7 / 


= 0 


and 


d<P 

dm 


d<t> y 
dm' J 



(2.23) 


with boundary conditions that are compatible with the choices of Dirichlet and Neumann boundary condi¬ 
tions. 

We specialize the Porous Medium model we are studying by choosing the second gradient part of the 
dimensionless potential energy, that is we assume 


,e',m, e) := -[ki(£') 2 + 2 k. 2 £'m' + k. 3 (m') 2 ]+ x ¥(m,£) 


(2.24) 


with ki ,&3 > 0, &2 £ M such that k i k$ — k 2 > 0. These parameters provide energy penalties for the for¬ 
mation of interfaces; they have the physical dimensions of squared lengths and, according with the above 
mentioned conditions, provide a well-grounded identification of the intrinsic characteristic lengths of the 


one-dimensional porous continuum. In this case, equations (2.231 become 


d'P S 

-r - (kie 1 + kim'Y = 0 and —- ikie' + ^ 111 ')' = -~-m. 

de dm pg f 


(2.25) 


We notice immediately that such a system of PDE has the form of the problem (P) introduced in Section 2.1 


provided the first gradient energy T* is a polynomial in the strain and in the fluid content. The main difference 
between the two system of equations lies in the fact that in the first of the two equations the derivative of 
the strain with respect to time is missing. This is due to the fact that, for simplicity and for coherence with 
the previous paper on which our sketch of derivation is based, we have not considered the dissipation forces 
acting on the solid components. If those forces would be taken into account, we would get a parabolic- 
parabolic system as the one in problem (P). 

An important application of the theory briefly recalled in this section is that to the study of phase transi¬ 
tions in porous media under consolidation, namely, when the system is acted upon by an external pressure. 
This issue will be discussed in Section [6] 


3. Proof of Theorem [3] 

In this Section, we prove Theorem [3] via a Leray-Schauder fixed point argument. Firstly, we study a 
regularised version of Problem (P) for which we prove the existence of a strong solution, see Theorem [8] 
The proof of Theorem[8]is divided in two main steps: In the first one, we introduce an auxiliary problem, 
depending on a parameter £ £ [0,1], for which a direct application of the theory of quasi-linear parabolic 
equations gives a unique classical solution. The second step is concerned with the definition of a nonlinear 
mapping that satisfies the hypothesis of the Leray Schauder argument, Theorem [9] Once Theorem [8] is 
proven, we exploit weak convergence methods to get the obtain the conclusion of the main Theorem [3] 


3.1. The regularized problem 
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Let us introduce the following mollified version of problem ( P ), namely: Find the pair (E.m) satisfying 


— +div(— k\Ve — fcV 5 /??) = fi(m,e) 

inQ, 

(3-1) 

^—bdiv(—— kfWm) =/ 2 (m,e) 
dt 

infl, 

(3-2) 

e(0) = £o, 

infl, 

(3-3) 

m( 0) = mo 

in £2, 

(3.4) 

£(h,t) = £ D (t) 

in S', 

(3-5) 

m(h,t) = m D (t) 

inS, 

(3-6) 


inS, 

(3-7) 


in S. 

(3.8) 


We refer to ( |3. 1 [ )— ( f378| ) as problem (P s ). In this section, we prove the existence of strong solutions to 
problem ( P s ). 

To define problem (P S ), we use the following definition of the mollified gradient of a function / [see 

e-g. (13)]: 


V 5 /:= V 


/ Js{x-y)f(y)dy , 

Jb& W J 


(3.9) 


where J§ denotes the standard mollifier defined for example in [ 1 ] and B s (x) is a ball centred in x E £2 with 
radius 8 > 0 chosen such that x + 8 E Q. A mollified function u enjoys of the following properties: 


Theorem 7. Let u be a function which is defined on M" and vanishes identically outside Q. If u G L p (£2), 
1 < p < °o, then J§*u E LP (fl). Also 

||^5*w||lp(£1) < ||w||lp(£1) an d lim \\J§ * u — «||iP(n) = 0. 


Also, for all f E L°°(Q) and 1 < p < it exists a constant c s > 0 such that 

l|V 5 /l|w ( n)<c 5 ||/|| L 2 (n) . (3.10) 

Note that as 8 —> 0, typically c s —> °o. 

Theorem 8. Assume Problem (P s ) has at least a strong solution 

(£g,m§) E £d + VL ' 1 (LI x S) x mo + Vf ' 1 (LI x S). 

Proof. The strong solution of the regularized problem (P 5 ) is obtained here by a direct application of the 
Leray-Schauder fixed point theorem, viz. 

Theorem 9 (Leray-Schauder Fixed Point Theorem). Let X be a Banach space and let PL be a completely 
continuous mapping ofX x [0,1] into X such that PP (y. 0) = 0 for all v E X. Suppose there exists a constant 
R > 0 such that 
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for all (v,C) F X X [0,1] satisfying v = 2F(y, £). 77m« t/ie mapping ofX into itslef given by 

^t(v) := ^(v,l) 

has a fixed point. 


A nice proof of the Leray-Schauder Theorem can be found e.g. in [ 121, Theorem 11.6; see also [18]. 

Solution to an auxiliary problem. For any given couple (e,m) with 

£ e Ce D +L 2 (S;V),m e £/n D +L 2 (S;V) 
and £ E [0,1], consider the initial boundary value problem 


d£ 

— +div(— k{We — kiV 5 m) = F\(e,m,e) 
at 

d m 

-^-+div(—fov e — kf\m) = F 2 (£,fh,m) 

e(0) = C^o, 
m(0) = £mo 
£(/i,f) = ££ D (f) 
m(h,t) = C,mn{t) 

d fd,,) = o 


in £2 x 5, 

in £2 x S, 

in £2, 
in £2, 
in S, 
in S, 

in S, 
in S. 


where we have set 






,(1) Jl) 


A(e,m,e) := £ + £ A^m^£ + £ £ *£ 


k2—0 k2—0 k\ —2 &2—0 

if |e|, |e| < Mi and |m|, |m| < M 2 and Fi(e,m,e) := 0 otherwise, and 


( 2 ) 


( 2 ) 


w (2) ^(2) 


F 2 (e,m,m) := £ A^e fe + £ A ( ki [e k 'm + £ £ A^g 

=0 k\=0 k2=2k\—0 


W oh^k 2 \ }n 


(3.12) 

(3.13) 

(3.14) 

(3.15) 

(3.16) 

(3.17) 

(3.18) 

(3.19) 


(3.20) 


(3.21) 


if |e|, |£| < M 1 and |m|, |m| < Mi, and Fi(£.m.m) := 0 otherwise. 

We note that F\ (e.m.e) = J] (e,m) and A(£■ m.m) = fjie.m). We split the proof of existence of solu¬ 
tions to the system ( |3.12[ )-( |3.19| ) into two steps: 

Step 1. The system of equations (|3.12[), (|3.14). (|3.16|) and (3.18) is a special case of the problem 


e t -gi(x,t)e xx -g 2 (x,t,e,e x ) =0 
e(x,0) = ££ 0 
e(h,t) = £e D (t) 
de 


in £2 x S, 
in £2, 
in 5, 


(3.22) 

(3.23) 

(3.24) 
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In our case, we have 


gi(x,t) =h, 
g 2 (x,t,e,e x ) = 


4‘> 




,r> 


4 1} 


k 2 =0 


k 2 =0 


ki=2k 2 =Q 


-k 2 div(V 5 m)+ £ A$m k2 + £ A^m* 2 e + £ £ A^m-e 


k 2 1, 


Under the assumptions H\ — H\ and trusting the classical theory of quasi-linear parabolic equations (see 
Theorem 7.4, Chapter V in [ 161), for any given 


(e , m) G L 2 (S; H 1 (ft)) x L 2 (5; H 1 (ft)) 


the problem (|3.12M3.14M3.16|) and (J3. 1 S)) admits the unique solution £ G V^"’ 1 (£1 X S). 

Step 2. The system of equations (|3. 13|),(|3.15J),(|3.17|) and (|3. 19) can be treated in an analogous way as in 

Step 1, in fact such system is a special case of the problem 

m t ~ £ 3 (x, t)m xx - gfixjpn. m x ) = 0 

in Q.X S, 

(3.26) 

m(x, 0 ) = £mo 

in £2, 

(3.27) 

m(h,t) = C m D lt) 

in 5, 

(3.28) 

dm,, 

*r ( ' 2i,)=0 

in 5. 

(3.29) 

In this case, we have 

g3(x,t) =k3, 

g 4 (x,t,m,m x ) = —k 2 div(V 5 e) — am 2 m — lab £ mm 

(2) 

- b 2 OC£ 2 m + £ A^e* 2 
k\ =0 1 



( 2 ) 


+ L A t-,i 

* 1=0 


^e kl m 


(2) (2) 


+ IL H A *ifa 

k 2 =2ki =0 


£ kl m kl l m. 


Based on the assumptions H 1 -H 4 and relying once more on the theory for quasi-linear parabolic equations 
(see Theorem 7.4, Chapter V in [16J), for any given 

(£,m)eL 2 (S-,H\D.))xL 2 {S-,H\a)) 


the problem ( |3.13| ),( f3TT5] ).( |3. 17] ) and ( 3.19 ) admits the unique solution m G V 2A (D. x S). 


Finally, we conclude that for any given couple (e,m) G L 2 (S',H l (Q)) x L 2 (S',H l (£1)) and 0 < C, < 1, we 
have (e,m) G V 2 2,1 (Qx S) x V 2 1 (Tl x S) as a solution to the problem ( |3. 1 2[ > ( [3 .1 9[ ). 

Definition 10. Denote by X := L 2 (S',H l (D.)) x L 2 (S\H { (Q)), fake arbitrary (e,m) G X, £ G [0,1] and take 
the couple (e,m) G V 2 x ,S) C X ns the solution to problem ( |3. 1 2[ )-( |3. 1 9[ ). We define the nonlinear 
mapping : X X [0,1] —> X by means of the equation 


{e,m)=&(e,m,Q. 


(3.30) 
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Basic a priori estimates. Now, we prove unifrom estimates for all (e,m) E X satisfying the equation 
(e,m) = for some C, E [0,1], In fact, note that if (e,m) E X is a fixed point of then 


(e,m) E V 9 2,1 (£l x 5). Define now 

Gi(e,m,e) := Fi(e,m,e) + d t £o, G 2 (£,m,rri )2 := F 2 (jt,m,m) +d t mo, 

and test ( |3.12[ ) by <p E V and ( |3.13| ) by y E V. 

We multiply ( |3.12[ ) by (p = e and integrate it over O x 5, to get 

1 d 

2 dt 


Consequently, we obtain 

\i! ns 


2 dt 

< 


J s + h J || Ve\\ 2 L2{cl) ds = 

~k 2 J s (^ 5 m,Ve) L 2^ ) ds + J^Gi(£,m,e),£) L 2^ ) ds. 
II L 2 (£l)d S + kl j 


(3.31) 


(3.32) 


(Gi(e,m,e),e) L 2 W ds + c^tf; I \\V s m\\ 2 L2{a) ds + ri I ||Ve||^ 2(n) ^ 


where we have applied the Young inequality on the right hand side of ( |3.31| ). It is easy to prove that there 
exist ci ,C 3 > 0 such that 


(Gi(£,m,£),£) L 2 ^ 


ds < Cl +C3 


; IIl 2 (£2)^ 5 - 


Thus, defining Cj := ci + c ri k 2 T\\'V 


5 112 


m 


iL~(n) 


, C 2 := (ki — r/), and then taking I] < k\, we obtain 


£ llL 2 (n)-H £/A " 12 


T ci J < Ci +C3 J | 

Applying Gronwall’s inequality, we get the following uniform estimate for £ 

II £ IIl 2 (S;V) < C, 

for a positive constant C independent of £ and S. 

We multiply now ( 3. 13| ) by ly = m and then integrate it over D x S, to get 

1 d I 

|| m 


£\\- L 2 {n) ds. 


(3.33) 


(3.34) 


(3.35) 


\\ 2 L2{n) ds + k 3 Js\\Vm\\ 2 L2{a) ds = 
kl /( v5g, V m )i 2 (n)^ + J(G2(e,m,m),m) L 2(fyds. 


(3.36) 


Consequently, 


\ d_ r 
2 dt Js 


sds < 


2 L2{il) ds + k 3 jf ||Vifi||| 2(0) < 

Js{G 2 (£,m,m)pn) L 2 ^ ) ds + c v kj Js\\V S £\\ L 2 {Q) ds + ri Jj\Vm\\ 2 L2(Q) ds, 


(3.37) 


where we have applied the Young inequality on the right hand side of ( 3.37 1 . Now it is easy to find c 4 > 0 
and C6 > 0 such that 


J (G 2 (£,m,m),m) L 2 ( a ) ds <c 4 + c 6 Jj\m\\ 2 L 2 ^ds. 


(3.38) 
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Take C 4 := C 4 + c T1 t^T\\ V 5 <j|li“(n)’ c 5 := (&3 — tj), and then taking 77 < k 3 , we obtain 


\ (ll B, lli2(o)-|l m (°)lli2 ( Q ) ) +C5^||Vm||2 2(n) ^<C4 + c 6 ^||m||^ ( ^. (3.39) 


By the Gronwall argument, we get the desired uniform estimate for m 


\\ m li 2 (S;V) — C) 

(3.40) 

for a positive constant C independent of £ and 8 . 

The following estimates are a direct consequence of (3.35 1 and (3.40 1 : 


ll e llt“(sx 2 (n)) — c > 

(3.41) 

m L°°(S;L 2 (I2)) < c ■ 

(3.42) 

Now, we observe that, by construction, actually f\ <G L°°(S’,L°°(Q.)). 

In particular f\ € L 2 (S\L?(Q.)), so 

we can show that 


L 2 (S;L 2 (ft)) ^ C. 

(3.43) 

The same property holds for m, i. e. 


\\^t m L 2 (S;L 2 (H)) ^ C- 

(3.44) 


Having established these basic estimates we are ready to complete the proof of Theorem [ 8 ] using the Leray- 
Schauder approach. 

Leray-Schauder fixed point argument. Take X := L 2 (S',V) x L?(S\V) and : X x [0,1] —> X defined by 
C ) = (e,m) (see Definition [To]), where (e,m) is the solution to the auxiliary problem < 13. 1 2[ >( f3 .1 9| ). 
First of all, let us prove that Sf : X x [0,1] — > X is continuous. To this aim we follow the spirit of [20]. 
For any sequence C„) e X x [0,1] such that 

(g„,m„,£„) -A (e,m,Q inX x [ 0 , 1 ], 

we denote by e„ the solution of the auxiliary problem 


(-^5 ^ 2(^5 £n,x) —0 

inD x S', 

(3.45) 

£ n (Xi 0 ) — l^n^O 

inG, 

(3.46) 

En(h,t) = C«e/j(0 

in 5, 

(3.47) 

r)f.. 


in 5, 

(3.48) 

ion of the auxiliary problem 

<? ( m„ -gi{x,t)m xx - g 4 (x,t,m,m x ) = 0 

inD x 5, 

(3.49) 

m„ (x, 0 ) = („m 0 

inD, 

(3.50) 

= C, n m D {t) 

in 5, 

(3.51) 


in 5, 

(3.52) 
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while (e,m) is the solution to the auxiliary problem ( |3. 12 1-( 3.19 1 . 

Subtracting the corresponding equations and testing with £„ — £ and m n — m, we obtain: 

(df£n — d t £,£ n —£)l 2 (Q.) + (div(— k\W £ n +ki'V£),£ n — £) 1,2 (ci) 

+ (div(-k 2 V 5 m n +k 2 V 5 m),£ n - £) i2(n) = (Fi(e„,m„,e) -Fi(e,m,£),£„- £) L 2 (n ), 


(3.53) 


and 


- d,m.m n - m) i 2 (n) + (div(-k 3 Vm„ + k 3 Vm), m n - m ) L2(n) 

+(div(-k 2 V 5 £„ +k 2 V 5 £),m„ -m) L 2 (a) = (. F 2 (£ n ,m„,m ) -F 2 (£,m,m),m n -m) L 2^y 


(3.54) 


2 ^ ||e »- £ lll 2 (£i)+^ll v ( £ «- £ )lll 2 (n) <’?ll v (e«-e)IIZ 2 ( a) 


Using Young’s inequality in ( |3.53| ), we get for rj > 0 

1 ^ 11 • 112 , ||T7/r. ^M |2 

+£ 77 k? 11 V 5 (m„ - m) \\ \ 2 + (/) (£„, m n , £) - iq (g, m,£),£„-£ ) i2 (n) . 
Now, it is straightforward to show that 

|A(e n ,w„,£)-Fi(e,w,£)| <ci|g n -g| +c 2 |m„-i 


- m 


(3.55) 


(3.56) 


The constants and c 2 can be computed here explicitly if needed. 

From ( 3.54 ) and proceeding in analogous way as for ( |3.53 ), we have for rj > 0 

\j t W m n ~m \\“ 2(n) + k 3 1 |V (m n -m)||^ 2(n) < T 711 V(m„ - m )\\ 2 l2(q) 

+ c n ^2 11 v<5 (e« - e) 11 L 2 (£2) + (^2 (e« , , m) -F 2 (e, m ,m), m n - m) i2 ( n) . 

It holds, as before, that 

-F 2 (g,m,/n)| < c 3 |g„ — g| +c 4 \m n —m\. 


(3.57) 


(3.58) 


Summing up ( |3.55[ ) and ( | 3.57 [ ) and using Young’s inequality and the property ( 3.10 ) in combination with 
Poincare’s inequality, we obtain 


1 d 

2 dt 


(l|e« - £ IIl 2 (£2) + IK -m||£2 (Q) ) + (*t — T7)||V(£„ - £)||\ 2{a) + (* 3 - r/)||V(m„ -m)||£ 2(Q) 

— ^t l|V(m„ — rh ) || L 2 ^ +C 2 II ^(^77 ~ IIZ 2 (fi) + c J? c 3p77 ^11 ( c 2 + C4) ||/n n — Hl 2 (£ 2 ) 

+r 7 / (£„ - £) 2 dx + ri / (m n — m) 2 dx. (3.59) 

We make use again of the Poincare inequality, so that we finally get 

(jl £ « — £ llz,2(^) + \\m n — m 1122(a)) II V(£„ — £) || i2 ^ +C 4 II V(m„ — m) ||£ 2 (q) 


1 <7 

2 dt 


< Ci||V(m„ -m)||£ 2(n) +C 2 ||V(g„ -g)||£ 2(fl) + TK 3 || £ » - g||£ 2(£2) + r](c 2 + c 4 )||m„ -m||^ 2 (a) . 


(3.60) 
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where C 3 := k\ — rj — t\Cc 3 , C 4 := £3 — rj — Crj (c 2 + cf) and C is the costant of the Poincare inequality. 

(3.61) 


From (3.60 1 , we obtain (in a compact form) 


X is compact. By the estimates ( |3.35[ ) and ( 3.4Q[ ), 


maps 


|| (£«) m ni C«) (^i m i 0 llx — c ( II i^ni^ni Cn) (fd^i C) llx ) ■ 

Thus the continuity of is proven. 

Let us now prove that W : X x [0,1] __ 

bounded sets from X x [0,1] into bounded sets of V ^' 1 (Q. x ,S). Since the embedding Vy ' x 5) into 
L 2 (S\V) is compact (compare ( |2.15[ ), cf. Aubin’s Lemma) then Sf : X x [0,1] —> X is also compact. 

Now, for £ = 0 we have 0) = (0,0) for all (e,m) G X. The esimates ( |3.35[ ) and ( |3.40[ ) imply that the 

solution of the problem (e,m) = C), for some £ G [0,1] is uniformly bounded in X. The existence of 

at least one fixed point to i.e. (£§,m§) G X with ^(Eg,mg, 1) = (Eg, mg), follows by the Leray-Schauder 
Theorem. Consequently, the equation (Eg,mg) = ^(£g,mg, 1) G V 2 ’ (£2 x .S) has a solution to Problem fP"), 
and so, Theorem[8]is therefore proven. 


‘IB: 


3.2. Passage to the limit 5—^0 

To complete the proof of the main result stated in Theorem [3] we pass now to the limit for 5—^0. In 
other words, we study the weak convergence of the solution (eg,mg) to problem (P 5 ) to the solution (e,m) 
to problem (P). 

We write down problem (P s ) in the form 


d ^ + div(-ki'V£g-k 2 V S mg) = fi (mg, Eg ) 

in £2 x S, 

(3.62) 

d +div( k 2 V S £g k 3 Vtng) - f 2 (mg,E S ) 

in £2 x S, 

(3.63) 

£g( 0)=£ 0 

in Cl, 

(3.64) 

mg( 0 ) =m 0 

in Cl, 

(3.65) 

£g(h,t) = S D (t) 

in S, 

(3.66) 

mg(h,t) =m D (t) 

in S, 

(3.67) 

£«■•*>=° 

in S, 

(3.68) 

t 5fe ')=° 

in S. 

(3.69) 


The next Lemma recapitulates the basic convergences we rely on. 

Lemma 11. Under the assumptions H 1 -H 4 , the following convergences hold up to subsequences, as 8 —> 0: 

(i) Eg —*■ £ in L 2 (S\V), mg —*• m in L 2 (S\V). 

(ii) d t £g —^ d t E in L 2 {S\L 2 {fl)), d t mg—*■ d t min L 2 (S’,L 2 (Q.)). 

(Hi) Eg £ in L 2 (S\L 2 (Q.)), mg m in L?(S\L 2 (Q.)). 

(iv) f\(Eg,mg) —> fi (£,m) a. e. in Q x S, fi(£g,mg) —> f 2 (£,m) a. e. in Qx S. 
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(v) V s £§ —¥ Ve in L 2 (S;L 2 (Cl)), V 5 m§ -^Vmin L 2 (S',L 2 (Cl)). 

(vi) £§ A £ in L°°(S;L 2 (Cl)), m§m in L°°(S;L 2 (£l)). 

(vii) Ifm§,£§ £ y 2 2+v,1 (n x S)for v > 0, then £ s ^ £ in L“(5;L“(«)), m 5 A m /„ L“(5;L“(ft)). 


Proof, (i) simply follows from the estimates ( ]3.35| ) and ( |3.40[ ), while (//') is a direct consequence of the 
estimates ( |3.43[ ) and ( |3.44[ ). To deal with (Hi), we make use of the Aubin’s compactness lemma (see Theorem 
[I}, particularly by choosing 

B 0 :=V, B = L 2 (kl)B i :=L 2 (Q.). 

Now, defining W (see ( |2.13[ )) as 

W := {(p £ L 2 (S-,V), d t (p £ L 2 (S'L 2 (Cl))} , 


we get 


W- 


■ L 2 (S',L 2 (Cl)). 


(iv) simply follows from (iii). (v) is a consequence of Theorem|7| (vi) follows from the estimates ( |3.4 1 [ ) and 
p.42[ ). To deal with (vii), we note that V 2 2+v ’ 1 L°°(Cl x S). In fact, we have H 1 +v (Cl) ^ L°°(Cl), V > 0, 


see [ 15] Theorem 5.7.8 p. 287, and hence, we deduce 


ll e 5l|L“(S;Z.“(n)) < C, 

11 11 z.-” (S;Z.“ (S2)) < C. 


From ( |3.70| > and ( ]3.71[ ) we obtain (vii). 


(3.70) 

(3.71) 

□ 


It is worth noting that the strong solution £ £ Ced + W' 1 (Cl x .S’), m £ Cm/) + Vf' (Cl x S) ensured by 
Theorem [8] is a solution of problem (P 5 ) and satisfies the identities 


>■ 2 , 1 , 


— / £g(x,0)(p(x,0 ) — / £§d,(pdxdt +k\ 

J£l JQ.XS 


+k 2 / V°m§V(pdxdt = , 

JClxS JClxS 


ClxS 

fi(£ S ,m s )(pdxdt, 


V £§V (pdxdt 


(3.72) 


— / mg(x,0)n/(x,0) — / mgdtXffdxdt +kj, 
Jo. JClxS 


+t <2 / V d £gVt ffdxdt = i 

JClxS JClxS 


Q.xS 

f2(£sd n s)¥dxdt 


VmgVydxdt 


for all test functions (p, y/ £ C(S;Cq(CI )) and c p(x,T ) = t jf(x,T) = 0 for all jc £ Cl. 


(3.73) 


The convergences (,)-(v) established in Lemma 11 are sufficient for taking the weak limit 8 —> 0 in (3.72 1 


and (3.731. Thus, we have 


— / £(x,0)(p(x,0) — / £d t (pdxdt + k{ / VEVcpdxdt 
JQ. JClxS JQ.xS 


+^2 / V°mV(pdxdt= / fi(£,m)(pdxdt , 
■JClxS JQ.xS 


(3.74) 
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— / m(x,0)yr(x,0) — / me), y/dxdt + k 2 / VmV y/dxdt 

JQ. JQ.xS JClxS 


+&2 


ISlxS 


V'V'V \j/dxdt = 


IClxS 


f 2 {£,m) y/dxdt 


(3.75) 


for all test functions (p and y/ previously chosen. 
Now, integrating back ( |3.74[ ) and ( | 3.75[ ), we obtain: 


{d t £, <P)l-(q.) + k\ (Ve, V 9 ) L 2 ( fl ) — —k 2 (Vm, V<p) L 2 ( n ) + (/i,<p)z 2 (n), (3.76) 

(<9 r m, t/r) L 2 ( fi )+ k 3 (Vm), Vt//) L 2 ( n ) = — k 2 {'Ve, Vy) t 2 (Q) + (/ 2 , V)l 2 (q.): (3.77) 

which is precisely our concept of weak solution to Problem (P), see Defintion[2j This completes the proof 
of the main result Theorem [3] 


4. Proof of Theorem |4| 

Let us suppose that problem (P) has two different solutions (£i,mi) and (£ 2 ,m 2 ) endowed with the 
same initial conditions. Define w e := £1 — £2 and w m := m\ — m 2 , then (w e ,w m ) satisfy, for all (p. ly G V, the 
following identities: 


{d t w e , (p) L 2 {a) + (div(-kj Vw e ), 9 ) L 2 (a) + (div(k 2 Vw m ), <p) L 2 (n) 

= (/i(ei,«i 1 )-/i(£2,m 2 ),(p) i 2( n ), (4.1) 

and 

V /r )z. 2 (ii) + (div(-fc 3 Vw m ), ^)l 2 (o) + (div(k 2 Vw e ), y) L 2 (n) 

= (/2(£i,"Ji) -h{ £ 2,m2), Y)l 2 {q.)- (4- 2 ) 

We choose (p := w e GV, y/ := w,„ £ V, then we obtain 

T^H We llL 2 (n) T || Vvre 11 i ,2 (o) = ~ k 2 j £i V w e Vw m dx + {}\ (£\.m\)- f\ (e 2 , m 2 )) w £ dx, 

(4.3) 


and hence, 


1 d 

2 dt 


w„ 


Il 2 (£ 2 ) ki 11 Vw 7 m 1 1 ^ 2 (q) = -k 2 J Vw e Vw m dx 

+ / (/ 2 (£ 1 ,/wi) — / 2 (£ 2 ,m 2 ))w m dx. 


(4.4) 


Applying the geometric mean-aritmetic mean inequality to the first term of the right hand side of both the 
above equations, we obtain 


2dt^ We ^ L2 {Q.) + ki\\Vw e \\ L 2( D: j < ^ ||Vw e || L 2( n ) + 2 l|Vw m || L2(n) 

+ / (/i(£i,«ii) - f\{£ 2 ,m 2 ))w e dx, 
Jci 


(4.5) 
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2 dt ll H ’ m ll i2 ( fi ) +^3||Vw m || i 2( Q ) < 2 ll^ w e|lz. 1 2 (n) + 2 ll Vw *lli?(A) 

+ / {f 2 (£i,mi)-f 2 (e 2 ,m 2 ))w m dx. (4.6) 

Ja 

We observe that, due to the particular structure of f\ and f 2 from ( |2.9| )-( [27TT[ |, there exist the constants 
A\, B\,A 2 , B 2 > 0 such that 


- fi(e 2 ,m 2 )\ < A^wel + Pi 


w„ 


|/ 2 (£i,mi)-/ 2 (£ 2 ,m 2 )| <A 2 \w e \ +B 2 \w m \. 
Now, we sum up (4A I and ( 4.6 1 and we use ( |4.7[ )-( |4.8| ) to get: 


1 d 


1 d 


2 ^H W e|| L 2 (n) + 2 ^IKIIl?(n) + (*1 -^ 2 )||Vw e ||£ 2{£2) + (* 3 -A 2 )||Vw m || L2 


2 dt 


(O) 


< +B\\w m \\ L 2^y 

Using assumption H 5 , we have 

2^H We ll L 2 (n) + 2d? ^ V "'^ 2 ( n ) - ^H We lli 2 (n) + H w '»llL 2 (a))> 

where A" = max {4 ,B\. Using now Gronwall’s inequality, we get for all t e S 

0 < \\ w e (0 IIZ 2 (t2) + H w '»(0llL 2 (a) — e (jl We (0)llL 2 (t2) + ll M, m( 0 )||i 2 , 

where w £ (0) = w e (0) = 0. Thus £1 = £ 2 , mi = m 2 a.e. in Cl for all t G S. 


(4.7) 

(4.8) 


(4.9) 


(4.10) 


(4.11) 


5. Negativity and boundedness of the strain and density: proof of Theorem [5] 

Due to the cross-diffusion-like structure the Problem (P) does not admit a weak maximum principle. 
Also, we were not able to find suitable test functions to obtain the boundedness or the negativity by means 
of an energy-like argument. In what follows, we rely on a regularity argument to ensure the boundedness 
property and use finite difference scheme to detect the negativity of the solution. 

The boundedness of the solution (£, m) to Problem (P) is obtained based on the additional regularity 
indicated in Lemma [IT] (vi) (together with 7/ ( ,). 

To prove the negativity of ( e,m ) we proceed as follows: Due to the regularity of the solution, in one- 
space dimension the functions £ and m are continuous. We then use the approach (and the corresponding 
notation) from [ 19] to construct using finite difference approximations a non-positive subsequence (e n .m n ) 
convergent to (£,m), situation which is valid provided the negativity of the functions f\. /? is guaranteed. 
We rewrite Problem (P) in the following way: 

1 It is worth noting that testing with (p = e+ e V and y/ = m + e V (or (p = (e — M\) + and \j/ = (m — AP)"* - to search for 

boundedness) does not work. This is mainly because we can not control the sign of terms like k 2 ^e~ Vm + or k 2 V^E~ Vm + . 
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d t m = (k 3 d x m + k 2 d x £) x +f 2 (£, m ) 

in £2 x 5, 

(5.1) 

5,£ = ki [d x e + kid x m) x + /i (£, m) 

in £2 x 5, 

(5.2) 

£(+ 0 ) = 6o(jc) 

in £ 2 , 

(5.3) 

m(x, 0 ) = mo(x) 

in £ 2 , 

(5.4) 

£(/i,0 =£/? 

in 5, 

(5.5) 

m{l\,t) = mo 

in S, 

(5.6) 

d x e(h,t) = 0 

in 5, 

(5.7) 

d x m{h,t) = 0 

in 5, 

(5.8) 


where we have f \, fi < 0 . 

Take now a positive integer N and let h := \/N. We introduce two types of grid points as 

Xj := ( i— ^)/z, i £ {0 ,,7V + 1} and jt; = ih, i € {—1,. ..,77+1} 
and set sub-intervals 

h := (xi-i,Xi), i e {1,.. .,N} and /, := (x ( -,x,-+i), i G {0,.. . ,7/+} . 
Furthermore, we set 


N 

I 

i=l 


*/, := £c,x, {c/}f = i Cl and W := C 


At 

L 

!=0 


where and % denote the characteristic functions of /, and I] n [ 0 , 1 ] respectively. 

Let {TnjJLj be a set of positive numbers and suppose that the n-th time step t n is determined by 


(5.9) 


(5.10) 


(5.11) 


/() — 0 , t n — tn — j + T n — Xfc, i £ { 1 ,..., ct} t m A T. 

k=i 

With such notation available, we find the unknown functions in the following form: 


for n = { 0 , 

Firstly, we introduce a difference scheme for equation (|5.2[>, viz. 


pii _ p n— t F » _1 - in n ~ l 

' ; =h — U i+[ +k 2 i ~ 2 '~ 1+ ' +/ 1 (ef- 1 ,/ W ?~ 1 ) 


tn 


h- 


h 2 

7 = 3,... ,7V and e„i = £i, £^+i = £w-t, 


(5.12) 


m h ■= L P 

i=l 

N 

a m(x,t „) 

(5.13) 

i=l 

N 

i £ v (x, t n ) 

(5.14) 

^Xi^F n = {hd x m 
i=0 

N 

+ kod x £) 

(5.15) 

4' 

i'=0 

~ e(x,L) 

(5.16) 


(5.17) 
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with £o = £d < 0 , £1 < 0 prescribed. 

Now we describe a difference scheme for ( |5.1[ ). We suppose that m? 1 is known from the time step t n -\. 
Then £^ can be calculated by ( |5. 1 7[ ). 


We compute b" 1 by 


b 1 ! = 


-i c (-1 


i E {0,... ,N + 1} . 


(5.18) 


Then we approximate the flux F" is by 


_ m n 

F, n = k 3-^- l -+k 2 bl i = — l 

n 

•N . 


Fq = 0 , F$ = 0 


(5.19) 

(5.20) 


Our proposed scheme reads then as follows: 

n— 1 


— m - ■ F n — F n , F " _1 — F "- 1 

' ' =0 1 . 1-1 +(1-0)-- ir -^+/ 2 (er 1 ,<" 1 ). 


/i 


h 




(5.21) 


with the boundary condition ( |5.20[ ), 9 E [0,1] and mo = m/> prescribed. 

Now we introduce the matrix representation of ( |5.22 1 and ( |5.20[ ). To this aim, setting A„ = t n /h 2 , we define 
the fV x N matrix H = [H^i] by 


-l, k = l = 1 , 

1 k = 1 , / = 2 , 

1 2<k<N — 1,1 = k—l, 

-2 2<k<N-l,l = k 

1 2 <k <N — 1,1 = k + 1, 

1 k = N,l = N- 1, 

-1 k = N,l = N, 

0 otherwise, 


and the A x A matrix B„ = by 


Bk.i ■= 


M2 • < 


—fci, k = l = 1, 

Z?2 k = 1, / = 2, 

0 2<k<N-l,l = k-l, 

-b k 2 <k <N —1,1 = k 
bk+1 2<k<N— 1,1 =k + l, 

b N k = N,l = N- 1, 
k = N,l = N, 


0 otherwise. 


(5.22) 


(5.23) 


Then ( |5.22[ ) and ( |5.20[ ) is reduced to 

(I - K 0H)m" - B„ 1 = (I + K (1 - 9 )H)m" - :1 + B„_, 1 + / 2 I1, 


« = { 1 , 2 ,, cr}, 


(5.24) 
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where m" = 7 (m”,.. -,m n N ), I denotes the identity matrix and 1 = l ( 1 ,..., 1 ). 

We state now some assumptions in order to prove the follwoing Theorem [12] which constructs teh sequence 
(, e n ,m n ) needed in the proof of Theorem [5] 

A\\ 2t„ 6 < h, 

A 2 : 2 t „(1 — 0) < h 1 , 

A 3 : b" -b^ x +b) -b}_ { <0. 

A \-At, are all technical assumptions. 

Theoreml2. Let n e {1,_a} and m n h 1 = Y,?=i £ i', 1 = E£Lo e f _1 & be given. Assume that 1 < 

0, e'iI 1 < 0 and m n h 1 , ej’ 1 are not identically constant. Assume A i—A 3 and define 

f h~ h 1 h~ 

p{e ’ h):=min \2(ree-y2el’ l( "' fa):= 2 S' 


Then if 


T n < rmn{p(e,h),l(h,k 2 )}, 


(5.25) 


t/ie scheme (5.22) and (5.20[) together with (|5. 1 7[) admits a unique solution 


N 


N 


Wrt) = CL e ?*i’L"tXi) 

i= 0 (=1 

which satisfies £ ( ” < 0, andm'j < 0/or all i G {0 ,..., A}, j G {1 ,..., A}. 

Proof. By assumption Ai, defining A := I — A„0H = [A*./] we can observe that the matrix A is irreducible 
and diagonally dominant. In fact the irreducibility is a consequence of 


Ak,k >0; 1 < k < N, 


(5.26) 


Aic,k—i <0; 2 <k <N, Ak t k +1 < 0, 1 < k < A — 1 . 

On the other hand, by assumption A\, we have 


(5.27) 


N 

> 1 

/=t 

for 2 < k < N — 1. In a similar way, we have 

N y N 1 

;=t z /=t z 

Thus A is diagonally dominant. From asssumption Ao, by a direct calculation it is possible to verify that 
every entry of (I — A„ (1 — 6 )H) is nonnegative. 

Now we write (|5.24|) in the following way: 


(I —A„0H)m" = (I + A„(l - 0)H)m"-' +B„_il + / 2 Il + B 1I l, 


(5.28) 
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and we exploit A 3 and fl5.25| ) to verify the negativity of the right hand side of ( |5.28[ ) and the negativity of £ ; ”. 
We procede by induction: 


is assumption A3. We suppose that 

br l -b’;:l+b';-bu<o 

holds and we prove the same inequality holds also for n. By ( |5.22[ ), we have 
b n_ b n_ i+b n+^ b n+l = 

bl-bt 1 + l (< +I -nil ~ ^ [*3« +1 - 2 m"+m"_ 1 )]) 

( T "^ [^K'+i 1 — 2 m ”- 1 +/<r 1 1 )]^ - (b'l-b’U) < 0 . 

Morover we know that /2 is negative. Thus, we have 

m" < 0 lor all I <i<N. 


(5.29) 


(5.30) 


(5.31) 

(5.32) 

(5.33) 


Take now n = 1 in ( |5.17[ ). Using both ( 5.33[ ) and the negativity of f\, we easily obtain £° < 0. Consider 
again ( |5.17| ) and suppose that £„_ 1 < 0, again by a direct calculation it holds £„ <0. □ 

Theorem [5] is now proven. 


6. Numerical study of steady states of strains and fluid densities for the consolidation problem 

As we already mentioned above, a very interesting application of the theory developed (recalled) in 
Section [T4| is the study of profile formation in porous media in a phase transition regime. We will consider 
a system exhibiting two phases differing in the strain £ and in the fluid content m. In this situation on a finite 
one-dimensional bar the system can show profiles, in £ and m, connecting one phase to the other. 

We consider the following expression for the total potential energy density in the perspective of describ¬ 
ing the transition between a fluid-poor and a fluid-rich phase 

^(m^e) := — m 2 (3m 2 — 8 b£m + 6b 2 £ 2 ) (6.1) 


where 


1,1 , 

T / b(/7?,£) := /?£ + -£ +-a(m — be) 


( 6 . 2 ) 


is the Biot potential energy density | 6 ], a > 0 is the ratio between the fluid and the solid rigidity, b > 0 is a 
coupling between the fluid and the solid component, p > 0 is the external pressure, and a > 0 is a material 
parameter responsible for the showing up of an additional equilibrium. 

In the papers [|8][9j we have studied the stationary version of the problem ( 2.25 [ ) corresponding to the 
potential energies (2.24) and ( |6.1[ ) to describe the possible occurance of an interface between two phases 
differing in fluid content. In fact this model is built in such a way to describe the existence of two states 
of equilibrium: the fluid-poor phase (£ s ,m s ) and the fluid-rich phase (£/,m/) corresponding to the two 
minima of the double-well potential energy l T in (|6.1|). Note that for a = 0.5, b = 1, a = 100 the pressure 
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ensuring the existence of two phases is p = 0.24221, while for the two phases we find e s = —0.1436, m s = 
-0.1436, e f = -0.1598, m f = -0.0427. 


The dissipative dynamics ( |2.25| ) of the porous medium model, starting from any initial state (go,mo), 
leads the system to a stationary state, which is the solution of the following problem: 


with 


(— k\e' —kjm')' 
{—kiE'—kpn)' 

e{h) = e D , 

m(l\) = m D , 




dm 

dx 


(h) = 0 . 


/t(m,e) 

/2(m,e) 


inD, 

inD, 


/t(m,e) := 

f 2 (m,e) := 


d x ¥ / 

dV . 
^ (m ' 


e) 

e) 


- bam 3 — ab 2 m 2 e — p — £ + abm — ab 2 e 
-am 3 +2abem 2 — b 2 ae 2 m — am+abe 


(6.3) 

(6.4) 

(6.5) 

( 6 . 6 ) 

(6.7) 

( 6 . 8 ) 


(6.9) 


where we recall ED- Note that this system of stationary equations has the same form of the stationary 


problem corresponding to our general problem (P) introduced in Section 2.1 


In our case of Dirichlet-Neumann boundary conditions, the stationary Problem ( |6.3[ )-( |6.8[ ) has not a 
unique solution. From the physical point of view, this property means that it is possible to observe different 
strain and fluid content stationary profiles with the same Dirichlet condition (edP^d) at one end. Below we 
discuss some graphs representing the solution of equations ( |6.3[ )-(6.8 1 in the interval = [0,1]. 

In figures 6.1 -6.2 the black solid lines correspond to the case k\ = &2 = = 10 3 , while the dashed 

lines correspond to the case k\ = kj = 10 3 and kj = 0.2 x 10 ' 3 , 0.8 x 10 3 . We recall that only in the 
second case the uniqueness of the solution to the time-dependent problem (P) is ensured (see Theorem 
[4]). We comment the physical features of the solution reffering to the fluid density profile m (bottom of the 
figures). The charachteristics of the strain profile can be discussed accordingly. 

The three different graphs for £ and m in figure |6.1| correspond to three different stationary solutions to 
problem ( ]6.3[ )—( [ 6 ~ 8 | ), for the same Dirichlet boundary value Ep = £, mp = m in /1 = 0, where £, m are values 
close to the fluid-poor phase (e s ,m s ) but slightly larger. In the top row the system is almost completely in 
the fluid-poor-phase, indeed in the interval [ 0 , 0 . 2 ] the profile quickly decay from m to m s and then it stays 
constant. Note that no interface between the fluid-poor and the fluid-rich phase. In the central row, after 
a quick transition from m to mj the system constantly stays in the fluid-rich phase. Even in this case no 
interface is seen. Finally in the bottom row the stationary profile is an interface between the two phases m s 
and mf. Indeed, the profile started at m first drops to m s and at a certain point quickly increases up to mp. 
The three different graphs for £ and m in figure |6.2| correspond to three different stationary solutions to 
problem ( |6.3[ )-( [N8] ), for the same Dirichlet boundary value Ed = £• = kn in l\ = 0, where £.m is the 

saddle point of the energy function <F. The essential features of the profiles are similar to those of figure [671] 
but the choice of the value m gives rise to small differences in the shape of the solutions. 

The solutions of the stationary problem ( |6.3[ )-( |6?8| ) are obtained numerically via the finite difference 
method powered with the Newton-Raphson algorithm. The use of different initial guess in the Newton- 
Raphson algorithm has allowed us to find numerically the different stationary solutions. In particular, for 
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Figure 6.1: Solutions e(jc) (left) and m(x) (right) of the stationary problem ( |6.3[ )—i [6~8] ) with the bound¬ 
ary conditions e(0) = £ = —0.141, m(0) = m = —0.13, d x e(x = 1) = 0, andd x m(x = 1) = 0 on the fi¬ 
nite interval [0,1], for p = 0.24, a = 0.5, b = l,a = 100,kj = 10~ 3 = 10 3 , k ,2 = 10 -3 (solid line), and 
k\ = 10 -3 = 10 3 , k .2 = 0.2 x 10 3 , 0.8 x 10~ 3 (dotted lines), starting by the following initial guesses (gray 
lines): costant fluid-poor phase (top), costant fluid-rich phase (middle), Dirichlet boundary conditions fix¬ 
ing the two phases at the ends of the sample (bottom). 
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Figure 6.2: Solutions e(x) (left) and m(x) (right) of the stationary problem ( |6.3| )-( [6?8| ) with the boundary 
conditions e(0) = £ = —0.1454, m( 0) = m = —0.0897, d x e(x = 1) = 0, andd x m(x = 1) = 0 on the finite 
interval [0,1], for p = 0.24, a = 0.5, b = l,a = 100, k\ = 10 -3 = 10 -3 , ki = 10 -3 (solid line), and k\ = 
10^ 3 = 10 -3 , k 2 = 0.2 x 10~ 3 , 0.8 x 10 3 (dotted lines), starting by the following initial guesses (gray lines): 
costant fluid-poor phase (top), costant fluid-rich phase (middle), Dirichlet boundary conditions fixing the 
two phases at the ends of the sample (bottom). 
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both figures 6 . 1 - 6.2 in the top row we used as initial guess a costant function equal to the fluid-poor- 
phase, while in the central row it has been used a costant function equal to the fluid-rich-phase. Fi¬ 
nally, in the bottom row we used as intial guess the solution to the same stationary problem with Dirichlet 
boundary conditions fixing the two phases at the ends of the sample, namely, (e( 0 ),m( 0 )) = (e s ,m s ) and 
(e(l),m(l)) = (£/,»?/) (see g). 

We now describe the adopted finite difference substitution rules. Let n be a positive integer number and 
let a = \/n be the space increment. We subdivide the space interval [0,1] into n small intervals of lenght a. 
Given a field h(x), for any i e {1,... ,n — 1}, we set 


h■('<*) ~ [*((*' + 1 )< 7 ) -*((*'- !)<*)] 


For the second space derivative we set 

h"(io,) « X^[h{(i+\)o)-2h{ic)+h{(i- l)cr)] 

for i G {1 ,... ,n— 1}. 


7. Conclusions 

We have studied the existence and the uniqueness of weak solutions to the problem (P) introduced in 
Section |2.1| We have stressed that the mathematical interest of this problem lies on the coupled cross¬ 


diffusion-like structure of the transport fluxes. The problem has, also, a remarkable physical application in 
the framework of the Porous Media theory (see the discussion Section |T4| ). 

It is worth noting that our mathematical approach is restricted to the one-space dimension case (as far as 
we are concerned with the passage to the limit 5—^0) and cannot be extended for higher space dimensions 
in a natural way. To make progress in this direction we hope to be able to employ the hidden variational 
structure of the problem [ fTOj Section 2.4], see also ( |2.18| ). Regardless the choice of space dimension, we 
find mathematically interesting the study of the t —> °° asymptotics in the case when multiple steady states 
are expected. Similar considerations can be made in the Cahn-Hilliard setup. 
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